Uniform description of bulk observables in the hydrokinetic model of A+A collisions 

at RHIC and LHC 



Iu.A. Karpenko 1 ' 2 and Yu.M. Sinyukov 1 
1 Bogolyubov Institute for Theoretical Physics, Metrolohichna str. 14b, 03680 Kiev, Ukraine and 
2 Frankfurt Institute for Advanced Studies, Ruth-Moufang-Str. 1, 60438 Frankfurt am Main, Germany 



m 
o 

(N 



(N 



K. Werner 

SUBATECH, University of Nantes - IN2P3/CNRS 



EMN, Nantes, France 



A simultaneous description of hadronic yields; pion, kaon, and proton spectra; elliptic flows; and 
femtoscopy scales in the hydrokinetic model of A+A collisions is presented at different centralities for 
the top BNL Relativistic Heavy Ion Collider (RHIC) and CERN Large Hadron Collider (LHC) 2.76- 
TeV energies. The initial conditions are based on the Glauber Monte-Carlo simulations. When going 
from RHIC to LHC energy in the model, the only parameters changed are the normalization of the 
initial entropy defined by the number of all charged particles in most central collisions, contribution 
to entropy from binary collisions and barionic chemical potential. The hydrokinetic model is used 
in its hybrid version (hHKM), which provides the correct match (at the isochronic hypersurface) of 
the decaying hadron matter evolution with hadronic ultrarelativistic quantum molecular dynamics 
cascade. The results are compared with the standard hybrid models where hydrodynamics and 
hadronic cascade are matching just at the non-space-like hypersurface of chemical freeze-out or on 
the isochronic hypersurface. The modification of the particle number ratios at LHC caused, in 
particular, by the particle annihilations at the afterburn stage is also analyzed. 

PACS numbers: 13.85. Hd, 25.75.Gz 
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I. INTRODUCTION 

Soon after the first Large Hadron Collider (LHC) 
heavy ion results were launched, it became evident that 
the hydrodynamic picture of the collision processes, con- 
firmed at Relativistic Heavy Ion Collider (RHIC),, is 
clearly seen also at much higher energy. This conclu- 
sion is based on two classes of observables. The first one 
is related to the azimuthal anisotropy of particle spectra 
expressed basically through its second harmonics, or v-i 
coefficients. The obtained LHC results for the transverse 
momentum dependence of 1*2 (pr) at a given centrality 
bin were found to be similar to the ones at RHIC, ex- 
cept for the higher momentum range at LHC pj. This 
is the evidence of the same hydrodynamic mechanism 
of the anisotropy formation as at RHIC. The second 
type of observables deals with the direct measurements 
of the space-time structure of nucleus-nucleus collisions 
by means of the correlation femtoscopy. The femtoscopic 
spatiotemporal scales of the systems are typically repre- 
sented in terms of the interferometry radii. The hydro- 
dynamic predictions Q for px behavior of the radii at 
the LHC energies were confirmed by the ALICE experi- 
ment 0. The most impressive hydrodynamic prediction 
[1, [f| , that the ratio of the two transverse interferome- 
try radii, out to side, will drop in the whole px- interval 
with increasing collision energy and reaches a value close 
to unity at the LHC, has been discovered experimentally 
1- 

Now the hydrodynamic-based approach to ultrarela- 
tivistic nucleus-nucleus collisions becomes the "standard 
model" for these processes. Note, however, that the de- 
tails of this approach did undergo essential modifications 



during the last decade and even now they are different 
in different models. The first attempts to describe suc- 
cessfully v 2 coefficients [|| at RHIC were based on the 
perfect hydrodynamics with an equation of state (EoS) 
corresponding to a first order phase transition. How- 
ever, a simultaneous description of V2, particle spectra 
and femtoscopy scales within the unified initial conditions 
collided with an invincible problem. Only in recent years 
have the main factors that make it possible to describe si- 
multaneously the spectra and femtoscopic scales become 
clear. They are |7HllT|: a relatively hard EoS because 
of a crossover transition (instead of the 1st order one) 
between quark-gluon and hadron phases, the presence 
of the prethermal anisotropic transverse flow developed 
to thermalization time, an "additional portion" of the 
transverse flow caused by the shear viscosity effect and 
fluctuations of the initial conditions. An account of these 
factors gives the possibility to describe well the pion and 
kaon spectra together with the femtoscopy RHIC data 
within a realistic freeze-out picture with a gradual decay 
of the fluid into observed particles [l2| ■ 

In fact, at the moment there is no real unique model of 
heavy ion collisions. Different hydrodynamic-based ap- 
proaches use different initial conditions as well as differ- 
ent final state treatments for the matter evolution, differ- 
ent EoS, etc. In this paper we make the next step towards 
converging to the "standard model" . We analyze a possi- 
bility for a uniform description of the soft observables in 
the hydrokinetic model of A+A collisions at RHIC and 
LHC, namely, a simultaneous description of the hadronic 
yield; pion, kaon, and proton transverse momentum spec- 
tra at different centralities; t>2 (pt) coefficients; and fem- 
toscopy scales in A+A collisions at RHIC and LHC. For 
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this aim we use the hydrokinetic model (HKM) [8l [L2l li~3l| 
in its hybrid version (hHKM) [l4| , which allows one to ap- 
ply hydrodynamics also at the late non-equilibrated stage 
of gradual system decay. It gives the possibility to switch 
correctly to the ultrarelativistic quantum molecular dy- 
namics (UrQMD) cascade at any spacelike hypersurface, 
in particular, at an isochronic one. We can compare the 
results also with a pure hybrid model where the hydro- 
dynamics and hadronic cascade are matching just at a 
non-space-like hypersurface of the chemical freeze-out or 
on the isochronic hypersurface. The goal of this paper 
is to fix the HKM by applying it to describe the known 
experimental data at RHIC and LHC and based on this 
to predict still unmeasured observables. The results pre- 
sented here are correlated with similar studies performed 
in other hydro based models, different from hHKM in 
various aspects [l~5l - [l7l ]. 

II. MODEL DESCRIPTION 
A. Initial conditions 

All the results presented in this work are calculated for 
the midrapidity region with the approximation of longi- 
tudinal boost-invariance. The hydrokinetic code [l2[ is 
developed now to simulate the full 3+ ID matter evolu- 
tion and it allows one to analyze the central, noncentral, 
and peripheral heavy ion collisions. The formation of 
the flow anisotropy depends on the initial eccentricity e 
which is defined by the centrality of collisions. Even at 
the maximal centrality the event-by-event fluctuations of 
the initial conditions (IC) lead to some nonzero effective 
eccentricity These fluctuations appear due to chaotic 
distribution of nuclcons (or local color charges) in col- 
liding nucleus and, as the result, the nucleon scattering 
positions fluctuate from event to event. In the hydro- 
dynamic picture such initial eccentricity of the (mean) 
energy density profile transforms into a flow anisotropy 
even at zero impact parameter. The effect persists for 
non-central heavy ion collisions, where it leads to some 
systematic increase of mean initial eccentricity. 

In this work, we employ the Monte-Carlo Glauber 
(MC-Glauber) ICs. 

1. The Glauber initial transverse profile 

The mean transverse density of nucleons in nuclei is 
defined as 

oo 

T ^= l u I 2 2 D I **!., (I) 

J exp[(Jx 2 L +x^- R A )/5] + l 

— CO 

where A is the atomic number; for Au+Au collisions one 
uses A = 197, then R A = 1.12A 1 / 3 -0.86A~ 1 ^ w 6.37fm, 
5 = 0.54 fm. Constant a is defined by the normalization 
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Table I: Cuts on the number of participating nucleons from 
MC-Glauber model, corresponding to different centrality 
classes for Au-Au and Pb-Pb collisions at top RHIC and 2.76 
TeV LHC energies, respectively. 



condition. In the so-called Glauber model the multiplic- 
ity (or entropy) produced in collisions of two nucleus is 
defined, roughly, by the number of participating nucleons 
contained in the overlapping distributions ([1]) for each nu- 
clei shifted by the value of the impact parameter. In the 
MC-Glauber approach, instead of utilization of the mean 
participant density ([T]), one starts the Monte-Carlo pro- 
cedure with Eq. ((TJ) as the distribution function of the 
random variable- nucleon number at point xt- The finite 
size of the nucleon can also be taken into account; this 
slightly modifies the average nucleon distribution com- 
pared to Eq. ([1]). The collision criterion for a pair of 
nucleons in each event is based on the value of nucleon- 
nucleon cross section gnn at the corresponding collision 
energy. The nucleon-nucleon collisions result in depo- 
sition of a certain amount of multiplicity (entropy) to 
different cells in a transverse plane. The contributions to 
multiplicity (entropy) from the "hard" elementary colli- 
sions and from the "soft" ones are different: The former 
are proportional to the number of binary collisions while 
the later are associated with the number of wounded nu- 
cleons, or participant number [lj§. The GLISSANDO 
code [lj| allows one to calculate the xt distribution of 
both numbers utilizing the MC-Glauber procedure. We 
use this generator and suppose that the initial entropy 
profile in the transverse plane at midrapidity is propor- 
tional to a linear combination of the density of wounded 
nucleons and that of binary collisions: 

/ \ ,l-5dN w dN hin , 

where C is the normalization constant, 5 = 0.14 for 
top RHIC energy and S = 0.08 for 2.76 TeV LHC en- 
ergy are fixed from the description of centrality depen- 
dence of multiplicity of all charged hadrons at midrapid- 
ity dN ch /dr] in a hydro+cascade model (2(1 [2l"j. 

The MC-Glauber model requires the (inelastic) 
nucleon-nucleon cross section a and the number of the 
nucleons N a ^ in the colliding nuclei as input. We fix 
this to (c = 42 mbarn, N &i h = 197) for Au-Au collisions 
for top RHIC energy, and (a = 64 mb, N a . h = 208) for 
the LHC case. Different centrality classes are defined via 
the corresponding cuts on the number of participants, as 
shown in Table |U The other choice is to introduce cuts 
on impact parameter instead (which is associated with 
a cut on average number of participants); however, we 
found that it does not lead to visually different profiles. 

Because we do not employ event- by- event hydrody- 
namics here, we average the profiles from a large-enough 
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200 GeV HH1C, in-plane 



200 QeV RHIC, out-of-plane 




Figure 1: (Color online) The initial energy density profile 
along transverse s/in-plane (left) and y/out-of- plane (right) 
directions in the classical Glauber model (black curve) versus 
the MC-Glauber model with variable-axes analysis (solid red 
curve) and fixed-axes analysis (dashed red curve). For the 
MC-Glauber case, 5 = 0.14 is used in Eq. (2). 



ensemble of events for each centrality class. The statis- 
tical fluctuations tilt in each event the principal axes of 
the ellipse of inertia and shift the center of mass rela- 
tive to the reaction-plane coordinate system. To account 
for this effect, we superpose the principal axes by ro- 
tation and rc-ccntering of each initial distribution and 
after that take averages over the ensemble of events (so- 
called variable geometry analysis, also implemented as an 
option in GLISSANDO code). This gives us nonzero ini- 
tial eccentricity as a consequence of fluctuations, even for 
simulations with zero impact parameter. The resulting 
averaged distribution is then associated with the initial 
entropy density for hydrodynamic evolution. 

In Fig. [1] we show the resulting energy density pro- 
files for 0-5% centrality collisions at the top RHIC en- 
ergy from the MC-Glauber model, compared to the pro- 
file from classical Glauber. Note that the second term in 
Eq. ([2]) corresponding to a contribution from binary col- 
lisions leads to a squeezed energy- density profile, which 
results in the increased radial flow in the system. At 
the same time, additional initial eccentricity (compared 
to the classical Glauber model) is gained mostly by per- 
forming the variable geometry analysis. 



2. Initial time and transverse flow 

In the hydrodynamic approach the initial entropy (or 
energy) density profile should be attributed to the ini- 
tial time Tj when (partial) thermalization is established. 
At this moment the system gains anisotropic prethermal 
transverse flow that is caused by the transverse finite- 
ness of the system 0, Q . Typically the flow is almost 
nonrclativistic, vt ~ yr, and for the Gaussian initial 
profile for freely expanding fields or partons, which are 
suddenly thermalized at the time r,, is proportional to 



the transverse coordinate r and inversely proportional to 
the homogeneity lengths squared 0] . If the homogeneity 
length is directed along an axis that is tilted by angular 
4> to the in-plane axis x, then one can write it in the form 



Vt = a 



(3) 



Naturally, the developing of the prethermal flows are 
different from that in hydrodynamic or free-streaming 
scenarios because of the different mechanisms of the flow 
generation at the prethermal and thermal stages. The 
initial velocity profile is also different from Eq. <j3j) be- 
cause of non- Gaussian initial energy distributions. The 
utilization of the method [22j to the thermalization pro- 
cesses in A+A collisions at RHIC and LHC, which are 
neither hydrodynamical nor suddenly thermalized free- 
streaming ones, is in progress now. However in this work 
we use a very rough approximation of ideal hydrody- 
namic evolution for the prethermal stage. This stage 
should start from the very initial time of the collision 
process, just after the c.m.s. energy in overlapping region 
of colliding nucleus exceeds their binding energy, say at 
ro = 0.1 fm/c; this time has to be used in a consistent 
approach. 

At initial time ro = 0.1 fm/c there is no collective 
transverse flow and our basic analysis is related just to 
this case. In addition, in Sec. IIII CI we demonstrate the 
possible corrections caused by the modifications of the 
final radial and elliptic hydrodynamic flows, which are 
qualitatively similar to a corresponding effect caused by 
the shear viscosity. 



B. Hydrodynamic evolution and equation of state 
in the chemically equilibrated zone 



The quark-gluon plasma and hadronic gas are sup- 
posed to be in complete local equilibrium above the chem- 
ical freeze-out temperature T c h with EoS at high temper- 
atures as in the lattice QCD. With the given ICs, the evo- 
lution of thermally and chemically equilibrated matter is 
described with the help of the ideal hydrodynamic ap- 
proximation. The latter is based on 3+1D (in the general 
case) numerical hydrodynamic code, described in Ref. 
[l2| . For this stage of evolution we use the latticeQCD- 
inspircd EoS of the quark-gluon phase [23j together with 
corrections for small, but nonzero baryon chemical po- 
tentials [TH, [H|, matched with chemically equilibrated 
hadron-resonance gas via crossover-type transition. The 
hadron-resonance gas consists of all (N = 329) well- 
established hadron states made of u,d,s quarks, including 
a meson [/o(600]. Quantitatively, the EoS table used is 
not visually different from widely used s95p EoS param- 
eterization [25], [26[, which is based on results from the 
hotQCD collaboration HI 111 . 
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C. Hydrokinetic and UrQMD cascade stages 

At the temperatures below T c h the system loses chem- 
ical and thermal equilibrium and gradually decays. In 
the hydrokinetic approach the dynamical decoupling is 
described by the particle escape probabilities in inhomo- 
geneous hydrodynamically expanding systems in a way 
consistent with the kinetic equations in the relaxation- 
time approximation for emission functions [l2[ • The 
method allows one, in principle, to take into account the 
back reaction of the particle emission on the system evo- 
lution, which is, however, a fairly complicated technical 
procedure that is equivalent to viscous hydrodynamics 
only at small deviations from local equilibrium, and we do 
not apply it. However, without accounting for the back 
reaction, a serious corrections of the emission function at 
large times are needed. It can be shown that the hydroki- 
netic approach without such corrections results in overes- 
timated effective temperature of proton spectra at RHIC 
energy, and not enough rise of interferometrv radii and 
volume from top RHIC to LHC energies Ml. To solve 
this problem we do not employ in this work the hydro- 
kinetic approach at large times until free-streaming of 
finally produced particles. Instead of this we switch over 
a hydrokinetic evolution of continuous medium to an evo- 
lution of particles within UrQMD hadronic cascade [2!| . 
The model which matches the HKM and UrQMD we call 
hybrid HKM (hHKM). 

1. The Cooper-Fry e switching over hydrodynamics to 
UrQMD cascade 

In the present analysis, we compare the two different 
approaches to match hadrodynamic and hadronic cas- 
cade stages. The first one is standard and corresponds 
to sudden transition from the hydrodynamic regime to 
the UrQMD cascade. In this approach, the distribution 
of i-th sort of hadrons at switching hypersurface is ex- 
pressed through hydrodynamic and thermodynamic vari- 
ables u(x), T(x), pi(x) by the Cooper-Frye formula [30| : 

PTdp T d(ppdy J 

• e^n^da^ (4) 

where n M is the hypersurface normal unit four-vector 
with a positively defined zero component n° > 0. The 
obtained distributions after the Monte Carlo procedure 
serve as an input for UrQMD cascade. It has been known 
for a long time that such a prescription leads to incon- 
sistencies [111, HIJ , if the switching hypersurface contains 



the non-space-like sectors, and so should be modified to 
exclude formally negative contributions to the particle 
number of the particles that move towards the fluid. A 
cut of the negative contributions 1 leads in its turn, to 
loss of causality and distorts hydrodynamic evolution. In 
addition, the particles interacting within UrQMD can in- 
teract also with fluid, and this opacity effect has to be 
taken into account, but usually not. So, despite its tech- 
nical simplicity, the Cooper-Frye prescription should be 
applied with caution. 

The method of direct matching of hydrodynamics and 
UrQMD cascade with Cooper-Frye prescription at the 
chemical freeze-out hypersurface we call traditionally a 
"hybrid" model. The chemical freeze-out isotherms typ- 
ically contain non-spacc-likc parts that can affect the re- 
sults for some observables, and we shall check how large 
the effect is, comparing the results with ones obtained 
at the hydrokinetic matching at an isochronic hyperesur- 
face. 

If the matching hypersurface that separates (pure) hy- 
drodynamics and UrQMD cascade is isochronic, we call 
this model a "hybrid-isochronic" one. In such a model 
the above-mentioned problems with non-space-like hy- 
persurface are absent. On the other hand, if one uses an 
isochronic hypersurface to match the cascade and hydro- 
dynamic stages, it can contradict to the obvious expec- 
tation that at the periphery of such a hypersurface any 
local-equilibrium input to the cascade stage is not cor- 
rect. The peripheral regions are spatially and temporally 
rather far from the freeze-out isotherm, have rather small 
temperatures and thus cannot be described hydrodynam- 
ically. We compare the results of the "hybrid-isochronic" 
model with the hHKM approach because the latter allows 
one to match UrQMD with a decaying hydrodynamic 
system, having nonequilibratcd particle distributions, at 
isochronic hypersurfaces and in this way to overcome all 
mentioned problems. 



2. Matching of the hydrokinetic stage and hadronic cascade 

To construct the matching between hydrokinctics and 
UrQMD one calculates the distribution function for each 
sort of hadron at the switching hypersurface. Using the 
technique of the integral form of the Boltzmann equation 
B Ell , one can express the UrQMD input as the collec- 
tion of all particles that propagate freely from the points 
of their last collision or just from the point of their cre- 
ation and do reach this hypersurface. Because the evolu- 
tion time is the Bjorkcn proper time t, the distribution 
is expressed naturally in hyperbolic coordinates. Then 
the input for UrQMD is constructed as: 
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\(T,x,p) = f l t - eq -(T ,x^(r),p)exp ^- J R^x^ (r),p)d,s j + 



(5) 



exp I - / [^(3,^^),^) + A(a,^ a) (r),p)]da 



r 



is the proper time, 



here x = {rr,0}, r = \A 2 - 
rriT = y m 2 + p\ is the transverse mass, 6 = n — y, i] 
is the space-time rapidity, and y is the particle rapidity; 
notation 



(r) 



--{sh~ 



t)} 
Pt 



(6) 



s/i(0)) ,f T - ^(Tche-y / s 2 + T 2 sh 2 0)} 



ttit 



(7) 



Eq. 

'(P,T) 



(O are as follows: 
is the collision rate of 



is used. 

The different terms in 

the i-th sort of hadrons with the rest of particles, 
G i ecay (X,x,p) is an income of particles into the phase- 
space point owing to resonance decays, and Di(X,x,p) = 
u c ™ h g Ti is the decay rate of a given resonance. To 
calculate the collision rates, we assume meson-meson, 
meson-baryon and baryon-baryon cross sections in a way 
similar to the UrQMD code |29| . 

The hydrokinctic approach allows one to build the 
(nonequilibrium) distribution function accounting for the 
above-mentioned effects on any hypersurface a sw that is 
situated behind the one corresponding to the hadroniza- 
tion process. As was discussed above in the consis- 
tent approach one should use the space-like hypersur- 
face to switch to the hadronic cascade. For the present 
analysis, in the case of hydrokinctic switching we use 
<7 SW : r = t sw = const hypersurface, which is completely 
space-like. However, other choices for switching space- 
like hypersurface are also possible. 

THe hypersurface of chemical freezeout, which is the 
matching hypersurface in pure "hybrid" scenario, is typ- 
ically non-spacc-likc, so the explicit paramctrization of 
the hypersurface t{tt) can be double-valued, as well as 
the paramctrization rrij, <f>). To escape multi- valued 
functions one may utilize these different paramctrizations 
in different sectors of the freeze-out hypersurface. In the 
case of the pure "hybrid" model it leads to the different 
representations for the distribution functions ^ for MC 
generation in different sectors of the hypersurface, 

d 6 Ni 

TTTTn TT - = &i ( r i Vi §i Ui PTi <Pp)t (8) 

drdrjdcpdydpT d<p p 

if the explicit dependence tt 

d 6 N, 

drx dndcfidydpx c 



Tt(t) is assumed, or 
Fi(r T ,?7,0,y,PT,0p), (9) 



if another form r = t{tt,4>) is used. For the isochronic 
switching procedure in the hydrokinctic approach, based 
on Eq. ([5]), only the form Eq. © is used. 

Once the distribution function F on the switching hy- 
persurface is known, we proceed with the standard (for 
present event generators) method of Monte Carlo event 
generation for input to UrQMD. First, one needs to cal- 
culate the maximum value of the distribution functions 
and the mean multiplicities of each sort of hadron in one 
event. Then, in each event we randomly generate the 
exact number of each sort of particle in a given rapid- 
ity interval according to Poisson distribution, such that 
the average multiplicity over a large set of events yields 
the mean multiplicity in the hydrokinctic approach at 
the chosen isochronic hypersurface. Then coordinates 
and momenta of each particle are generated randomly 
according to the distributions ((4]) or ([5]). To simplify 
the procedure of generation, the accept-reject algorithm 
is used, which means the following, e.g. in the case of 
distribution ©: We generate the random set of coordi- 
nates x r = {t, 77, (p} and momenta p r = {y,pT,4>p} an d 
a random variable a in interval a G [0 . . . -F ma x] , where 
^max is the maximal value of distribution ([5]). Then, 
we compare the value of distribution in the given point 
F r = Fi(x r ,p r ). If a < F r , then the particle coordinate 
and momentum set are accepted; otherwise it is rejected, 
and the new random iteration is made. In this way, we 
generate in each event the coordinates and momenta of 
all particles with \rj s \ < 2 and \y\ < 2 intervals, which are 
wide enough to study particle production at midrapidity. 

The generated set of particles is then used as input for 
UrQMD code (29|, which computes further elastic and 
inelastic collisions and decays of unstable hadrons. Fi- 
nally, the output of UrQMD - which is a set of particle 
positions and momenta - is analyzed to obtain physical 
results of the model. 



III. THE OBSERVABLES IN THE MODEL 

The azimuthally averaged transverse momentum spec- 
tra, 

dN 



2irpTdpTdy 



is calculated in the direct way from a generated event en- 
semble. When comparing to experimental data, we retain 
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or neglect the feed-down from weak decays in accordance 
with experimental procedure. 

In the described Monte Carlo procedure of event gen- 
eration, the elliptic flow coefficients are calculated in the 
model in a "standard" way, according to: 

2 2 

v 2 =< cos2(0 p - *rp) >=< Px ~ Py 2 >, (10) 

where the brackets correspond to the mean value over the 
particles from a given event ensemble (in a given pr bin, 
for pr-differential 1)2)- The procedure of reaction plane 
constructing is described in Sec. IIA1. 

For femtoscopy analysis, we first calculate the three- 
dimensional two-pion correlation function as a function 
of relative momentum q = p\ — p%, for each kx bin, where 
k,T = ij>\ +P2V2 is the average momentum of the pion 
pair. Following the experimental cuts (which are some- 
what different for STAR, PHENIX and ALICE collabo- 
rations), we consider pions in central pseudorapidity re- 
gion \q\ < 0.5. To calculate the correlation function, we 
use the same technique as in our previous studies in the 
FASTMC event generator [34|. Namely, the correlation 
function for bosons in the Monte Carlo procedure is cal- 
culated according to 

E 5A(q-Pi+Pj)(l + cos((pj -Pi)(xj - Xi)) 

(11) 

where 6a(x) = 1 if \x\ < Ap/2 and otherwise, with Ap 
being the bin size in histograms. We decompose the rela- 
tive momentum q into (ffout, <7side, QWg) projections, and 
preform analysis in the longitudinal center of mass sys- 
tem (LCMS), where the mean longitudinal momentum of 
the pair vanishes. 

Evaluation according to Eq. (fTTj) is done with the help 
of 3D histograms, implemented in ROOT library classes 
[35T ]. and two separate histograms arc used to calculate 
the numerator and denominator of (|lip . which are di- 
vided to get the correlation function. We fit the resulting 
correlation function with the Bertsch-Pratt parameteri- 
zation: 

C{q) = 1 + A • exp(-i^ ut g£ ut - R 2 idc q 2 idQ - R? ong qi ong ) 

Owing to longitudinal boost invariance and approxi- 
mate azimuthal symmetry for the most central collisions 
which we consider for the present HBT studies, the cross- 
terms R 2 s ,R 2 sl and are neglected. 

A. The results for RHIC 

The model parameters for MC-Glauber ICs are defined 
as follows. We consider the three switching procedures 
between hydrodynamics and the hadronic cascade. In 
the "hybrid" model the chemical freeze-out isotherm with 



temperature T c h (defined below) serves as the switching 
hypersurface. In "hybrid - isochronic" model the switch 
of hydrodynamics to the UrQMD cascade is implemented 
at the isochronic hypersurface t sw = const, and t sw is de- 
fined from the requirement that the isochronic hypersur- 
face touches the chemical freeze-out isotherm at r = 0: 
T(r sw ,r = 0) = T c h- In the hydrokinetic procedure the 
switching hypersurface is the same isochronic one but 
at time r < t sw the evolution is described by the hy- 
drokinetic but not not pure hydrodynamics. The MC- 
Glauber initial conditions are attributed to the starting 
time To = 0.1 fm/c in all the scenarios. The basic results 
are related to zero initial transverse flow. 

The normalization constant C for the initial entropy 
density in Eq. (|2|) is defined by the all charged particle 
multiplicity in most central collisions c = — 5% and is 
fixed for the given collision energy. The same concerns 
the parameter S that defines the binary collision contri- 
bution according to Eq. (0). For RHIC S =0.14. We 
also fix the parameters of chemical freezeout: T c h = 165 
MeV, [i B = 28 MeV, fj, s = 7 MeV, /ze = -1 MeV, ac- 
cording to the analysis of particle number ratios from 
STAR in thermal model [3(1 [37| , and include contribu- 
tions from weak decays to the proton spectrum to make 
correspondence with STAR procedure. Some of the re- 
sulting particle number ratios are shown in Table [TT1 Note 
that in our hadron table we also include additional res- 
onance states, e.g., Jo (600), /o(980), a (980) and high 
mass resonances with m > 2 GeV, followed by the recent 
compilation from Particle Data Group [H|. /o(600) is 
observed as a broad resonance structure, with very little 
knowledge about its decay channels and branching ra- 
tios. However, these resonances contribute to both final 
pion and (less) proton yields [36j; thus, modifying par- 
ticle number ratios, in particular, improves (anti)proton 
yields at RHIC. 



1. The yields, spectra, V2 and femtoscopy 

We focus first on the description of the top RHIC 
Au+Au data in hHKM model, which is shown in Fig. 
® 

It is known that even for the most central Au+Au col- 
lisions STAR and PHENIX proton multiplicities differ by 
a factor of w 2, pion multiplicities by more than 15%, and 
these discrepancies rise from central to peripheral central- 
ity classes, see Fig. [5] Within our model this difference 
cannot be fully reproduced by switching on/off weak con- 
tributions only, e.g. to reach lower pion multiplicity at 
PHENIX one also has to decrease initial energy density. 
We choose ICs in the model to reproduce the STAR mul- 
tiplicities, but compare the results with the spectra and 
HBT radii measured by both collaborations. Because 
transverse momentum spectra of identified hadrons are 
measured in a wide px range only by PHENIX, we mul- 
tiply the PHENIX proton spectrum by the factor 2 to 
better compare its shape with the STAR points and the 



7 



60-70% hybrid - 

60-70% hHKM - 
70-80% hybrid 

80-92% hybrid - 



0-5% hybrid 
0-5% hHKM 
0-5% isochrone 

— 20-30% hybrid 

- - 20-30% hHKM 

40-50% hybrid 
40-50% hHKM 




1.5 2 
P T [GeV] 



protons | 



60-70% PHENIX x2 
70-80% STAR 
70-80% PHENIX x2 
80-92% PHENIX x2 




1.5 
P T [GeV] 



10 £ 



" 10 
> 



* 1 - 



-10" 



* 10' 



lo- 



0-5% STAR 
0-5% PHENIX 
20-30% STAR 
20-30% PHENIX 
40-50% STAR 
40-50% PHENIX 
60-70% STAR 




* 60-70% PHENIX 

* 70-80% STAR 

a 70-80% PHENIX 
80-92% PHENIX 

I ' 11 ' 11 1 I ' 1 ' I 



0.4 0.6 0.8 1 1.2 1.4 
P T [GeV] 



V 2<P T ) a" charged 




0.2 0.4 0.6 0.8 1 1.2 
P T [GeV] 



7 - 



■0.1 



hybrid, c=0-5% 

hHKM, c=0-5% 

isochrone, c=0-5% 
hybrid, c=20-30% 




0.2 



0.3 



0.4 0.5 0.6 
P T [GeV] 



_l I I 1_J I I 



0.7 



0.8 



0.9 



PHENIX, Tt'Ti-, c=0-30% 
STAR, uriC, C=0-5% 
STAR, Tt'JI", C=20-30% 
STAR, Tt'JI", C=30-50% 




0.4 0.5 0.6 
P T [GeV] 




P T [GeV] p T [GeV] 

Figure 2: (Color online) ir~ ,K~ and proton spectra and elliptic flow of all charged hadrons at midrapidity for 200A GeV RHIC 
energy and different centrality classes (top) and HBT radii of tt~ pairs for most central events, calculated in different models 
and compared to experimental data from the STAR [H, [3, [|f| and PHENIX [H, collaborations. Note that proton 
PHENIX spectra is multiplied by the factor 2 (see text). Dashed lines correspond to the hydrokinetic procedure of matching 
(hHKM), while solid lines stand for the "hybrid" model case. The results for the "hybrid-isochronic" model are presented 
for central c = 0-5% events with a gray solid line. The ICs used: to = 0.1 fm/c with zero initial transverse flow and the 
MC- Glauber profiles for initial entropy density. 
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hHKM calculations. From Fig. [2] it is seen that the kaon 
multiplicity is over-predicted in hHKM for non-central 
events if compared to the PHENIX data, which could be 
a sign of the incomplete equilibration of strangeness in 
peripheral collisions 4l|. However, this effect is not seen 



when one compares results with the STAR kaon data. 





K- 


N K - 


N p 


Np 


STAR 


327T33 


49.5T7.4 


34.7T6.2 


26.7T4.0 


hHKM 


330 


47.3 


29.9 


20.5 



Table II: Identified hadron yields measured by STAR collab- 
oration [3^| for most central (0-5%) AuAu collisions at top 
RHIC energy, compared to hHKM calculations. 

In the middle row of Fig. [2]we compare the elliptic flow 
coefficients in full hHKM and "hybrid" scenarios with 
the STAR Q and PHENIX results for all charged 
hadrons, obtained with event-plane method. In addition, 
in Fig. [3] we present elliptic flow of charged pions, kaons 
and (protons+antiprotons), calculated with the same pa- 
rameters as for Fig. [3J and compared to i>2{2} for charged 
pions, kaons, and anti-protons from STAR (using a two- 
particle cumulant method). This calculation is done only 
for the centrality c=20-30%, because a considerably big- 
ger set of events must be used for this observable. One 
can see some overestimate of the data for ^-coefficients, 
it is because we use ideal hydrodynamics, but not a vis- 
cous one. 

The comparison of intcrferomctry radii, calculated in 
the model with the experimental data is shown in the 
bottom of Fig. H Note that PHENIX presented its re- 
sults for the 0-30% centrality bin, which corresponds to 
a smaller average multiplicity than the 0-5% STAR bin, 
thus PHENIX radii lie slightly below ones calculated by 
STAR; in our model we observe the same tendency with 
average ICs for 0-30% centrality. 

From Figs. [5] and one can conclude that both 
the hHKM and the "hybrid" models describe the data 
quite satisfactory, except 1)2 for very noncentral collisions 
(where the viscosity should play an important role and 
reduces elliptic flow) and HBT radii for 40-50% central- 
ity, which are clearly underestimated. The shear vis- 
cosity should also reduce the R ou t/Rside ratio because 
it enhances the transverse flow. As for the "hybrid 
- isochrone" model, it fails to describe the shape of 
pion, kaon and proton transverse spectra, V2, and long-, 
side- and out- interfcreometry radii. The main differ- 
ence between the first two (hHM, hybrid) and "hybrid - 
isochrone" scenarios is that the first two matching proce- 
dures do not use the local equilibrium particle distribu- 
tion functions as input for UrQMD cascade at the space- 
time regions where the system should be far from equi- 
librium. The peripheral regions at isochronic hypersur- 
face are spatially and temporally distant from the freeze- 
out isotherm and have rather small temperatures and 
in this transition area the finite and rapidly expanding 
system cannot be described hydrodynamically: The free- 



streaming regime of particle propagation already starts 
there. However, such a system can be described with the 
hydrokinetic approach. In what follows we consider only 
hHKM and "hybrid" models. 
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Figure 3: (Color online) Elliptic flow of tt ,K and 
(anti)protons at top RHIC energy and centrality c=20-30%, 
calculated in the hHKM and "hybrid" models with the same 
parameters as in Fig. [2] and compared to the STAR data [3| 



B. The results for LHC 

We now move to LHC energy to inspect now well the 
hHKM and "hybrid" models fit to 2.76 TeV LHC Pb+Pb 
collisions. First we have to account for the correspond- 
ing increase of the initial entropy production at the LHC 
energy by changing the normalization constant C in Eq. 
@ of the initial ener gy density profile, the parameter S 
is chosen to be 0.08 [21| for LHC energy. With bigger 
initial energy density, hydrodynamic evolution evidently 
leads to larger effective volume at the chemical freeze- 
out hypersurface, as well as to bigger transverse flow; see 
Fig. [4] The second modification is related to the chemi- 
cal composition at freeze-out: According to approximate 
particlc-antiparticlc symmetry at ^fs — 2.76 TeV energy, 
confirmed by ALICE data [4fJ, we set all chemical po- 
tentials to zero: [Lb = Hq = ^ s = 0. 



1. Yields and particle-number ratios 

Let us start with the particle-number ratios obtained 
in the hHKM model for LHC energy. In the simula- 
tions we also observe antiparticlc/particle symmetry in 
tt~/it + , K~/K + , p/p ratios, which are all close to 1. 
From the spectra plots in Fig. [5] one can conclude that 
(anti)proton multiplicity at midrapidity in the hHKM 
model is only slightly overestimated. Indeed, the non- 
trivial particle number ratios at midrapidity are shown in 
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v T [c] 



Figure 4: (Color online) Chemical freezeout hypersurface pro- 
files and transverse flows at the chemical freezeout position in 
the hHKM model. Profiles are made in the in-plane direction 
(solid lines) and the out-of-plane direction (dashed lines) for 
top RHIC (black) and 2.76 TeV LHC (red) A+A collisions. 
Initial conditions are taken from the MC-Glauber model. 

Table Hill To understand which factors contribute to suc- 
cessful description of particle- number ratios at LHC (and 
in particular the p/ir ratio), let us calculate hadron yields 
in different scenarios of evolution at the post-chemically 
equilibrated phase, while keeping the same ICs and chem- 
ical composition at chemical freeze-out. In this subanal- 
ysis we look at the most central collisions, where rescat- 
tering effects in cascade - via UrQMD code - should 
be most prominent. From Table IIIII one can see that 
both pion and proton yields are minimal for the "ther- 
mal model" scenario, where only resonance decays are en- 
abled. Involvement of UrQMD to calculate both clastic 
and inelastic scatterings (except for baryon-antibaryon 
annihilation, turned off with CT0ption(19)=l) increases 
somewhat both pion and lambda yields. Finally, switch- 
ing on £?i?-annihilation suppresses baryon yields, and at 
the same time increases pion yield, thus lowering the p/ir 
ratio to the value 0.052, which is quite close to the one 
measured by ALICE [40fl , Thus, one can conclude that 
annihilation processes in UrQMD are essential for suc- 



N-n N K N p Na p/tt K/ty A/tt 

full 775 123 40.5 20 0.052 0.158 0.026 

full-SB 716 114 50.5 24 0.072 0.159 0.034 

thermal 679 127 54 20.3 0.08 0.188 0.03 



Table III: Particle multiplicities and particle number ratios, 
calculated within the hHKM model for most central (0-5%) 
PbPb collisions with ^/s = 2.76 TeV in difieren scenarios of 
particle production: full scenario (hydro+UrQMD), full-SB 
(baryon-antibaryon annihilator switched off in UrQMD), and 
thermal model (kinetic phase with resonance decays only). 



ccssful reproduction of the p/ir ratio at the LHC energy. 
When going to nonccntral collisions, the p/ir ratio slightly 
increases up to 0.058 for 20-30% centrality (consistently 
with ALICE data), which should be the result of fewer 
inelastic processes owing to decreased effective volume at 
the hadronization hypersurface. Because charged hadron 
density at midrapidity for 20-30% central PbPb collisions 
at 2.76 TeV is close to the one for most central collisions 
at top RHIC energy, one can also conclude that inelastic 
processes in the cascade play a similar role also at RHIC, 
modifying particle number ratios and, in particular, de- 
creasing the p/ir ratio. It is worth noting that the value 
of the effect depends on the dynamics of the fireball that 
defines a duration of the afterburn stage and so can differ 
at RHIC and LHC energies. 



2. Spectra, elliptic flow and femtoscopy 

In Fig. [S]we show the comparison model/experiment 
for 2.76 TeV LHC energy. The experimental data are 
related to elliptic flow coefficients from the four-particle 
cumulant method, t>2{4} [H[ and interferometry radii to 
the most central events measured by the ALICE Collab- 
oration Q. In addition, as for the RHIC case, in Fig. 
IH1 we compare the elliptic flow v 2 {2} of charged pions, 
kaons and antiprotons for the centrality 10-20%, mea- 
sured by ALICE (49|. On the last plot, the experimental 
points arc somewhat above the hHKM calculations, and 
the probable source of the difference is the two-particle 
cumulant method, which give systematically bigger ellip- 
tic flow than the four-particle cumulant method used for 
the analysis of the elliptic flow of all charged hadrons. 
To compare with the experimental antiproton spectrum 
from ALICE, we exclude all feed-downs from weak de- 
cays, except for S + . As a result, we observe some dif- 
ference compared to ALICE for the mass dependence of 
the effective temperature (inverse slope) of the transverse 
momentum spectrum: the resulting antiproton and kaon 
spectra are too flat in the model, while the pion spec- 
trum is reproduced much better in a wide pr region. 
The reason for such a mismatch could be that the im- 
itation of the viscous effects by the initial transverse flow 
is not so effective at LHC than at RHIC because of a 
more protracted (viscous) hydrodynamics stage. Con- 
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P T [GeV] p T [GeV] 

Figure 5: (Color online) n~ ,K~ and proton spectra and elliptic flow of all charged hadrons at mid-rapidity for 2.76 TeV LHC 
energy and different centrality classes (top) and HBT radii of tt~ pairs for most central events, calculated in the hHKM model 
and compared to experimental data from ALICE 0, l4(il . |48| . Dashed lines correspond to hydrokinetic procedure of matching, 
while solid lines stand for the "hybrid" model case. Corresponding HBT radii for top RHIC energy are shown for comparison 
purposes. Note that calculations for c = 10 — 20 % are not provided, but the i^-coefficient for c = — 5 % is presented instead 
as prediction. 
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Figure 6: (Color online) Elliptic flow of n ,K and 
(anti)protons at top 2.76 TeV LHC energy and centrality 
c=30-40%, calculated in the hHKM model with the same pa- 
rameters as in Fig. [5] and compared to ALICE data [49l ] 



corning interferometry results, in hHKM we get system- 
atically lower i? s ide and i?i on g radii than the ALICE data 
for larger px', however, they are within the experimental 
error bars almost in all px regions. Note that the rise of 
interferometry volume observed by the ALICE collabo- 
ration is well reproduced in hHKM (see also Ref. [HI). 
Because we keep untouched the main model parameters 
when passing from RHIC to LHC energies, except for a 
general normalization for increased dN/dy, contribution 
from binary collisions, and the baryonic chemical poten- 
tials at freeze-out, one can conclude that the soft physics 
at RHIC and LHC is similar. 



C. Possible corrections to the basic results 

The above results are presented for perfect fluid dy- 
namics of quark-gluon plasma and noncquilibrium evolu- 
tion of hadronic gas after the chemical freeze-out. As it is 
known, quantum theory gives the limitation on the shear 
viscosity coefficient r\, more exactly, on the minimal ra- 
tio 2 = where s is entropy density. The inclusion of 
the shear viscosity in the boost-invariant hydrodynamic 
models affects the final hydrodynamic flow in the fol- 
lowing way: it increases the radial and decreases elliptic 
flows. 

In this work we do not utilize viscous hydrodynamics. 
However in this Section we try to demonstrate the influ- 
ence of the above-mentioned effects on observed trans- 
verse spectra, elliptic flows, and HBT. Of course, the 
conservation laws have to be taken into account to see 
reliable correlation in modifications of different kinds of 
observables. For this aim we use our model with modi- 
fied ICs that affect the final flow in the above-mentioned 
way. We found that the artificial addition of a small 
initial transverse velocity with profile ([3]) can serve for 



this aim. In this case the pressure gradient will drop 
faster in all directions and this reduces an efficiency 
of the transformation of the anisotropy of initial pres- 
sure gradient into the final flow anisotropy Q. For- 
mally, when v — > v + v rac i, the anisotropy parameter 
e v = (v x - v y ) I (v x + V y ) drops. 

From the comparison with the LHC experimental data 
in FigjT] one can see some improvement of the behav- 
ior for all the observables, when small initial flow is 
included with the average value being (ux,ini}=0.014 - 
0.018. Namely, the effective temperature of the single- 
particle transverse spectra at large pr increases, while 
the ^-coefficient decreases and i?out/-Rsidc ratio becomes 
lower. 

These "viscous-like" corrections bring the hope that 
the good description of the RHIC and LHC experiments 
can be reached within the HKM based on viscous hydro- 
dynamics. This will be the next step in our study. 



IV. CONCLUSIONS 

We present a consistent description of the transverse 
momentum spectra for the most abundant hadrons (n, 
K, p), elliptic flow coefficients, and interferometry radii 
for central and noncentral collisions at the top RHIC and 
LHC energies in the hybrid hydrokinctic model (hHKM). 
The latter provides the correct matching of decaying 
hadron matter evolution with hadronic UrQMD cascade 
at isochronic (or space-like in general) hypersurface. The 
results are compared with ones obtained at the different 
matching procedures in hybrid models with the same ini- 
tial conditions. It is found that the matching of hydro- 
dynamics with cascade at the chemical freeze-out hyper- 
surface gives the close results to hHKM, while the results 
are essentially different from both of these models and 
from experimental data when the matching of thermally 
and chemically equilibrated hydrodynamic systems with 
UrQMD happens at the isochronic hypersurface. 

The initial conditions allowing simultaneous descrip- 
tion of the soft observables in ultrarclativistic A+A col- 
lisions are associated with the MC-Glauber initial en- 
tropy density distribution with mixed contributions from 
wounded nucleon density and density of binary scatter- 
ings. One can see the satisfactory description of the 
RHIC and LHC data within perfect hydrodynamics for 
the QGP phase, except for some overestimation of v-2 
coefficients and R ou t/Rside ratios caused, probably, by 
neglect of the shear viscosity. As it is known the latter 
increases the final transverse flow and damps the ellip- 
tic flows. Qualitatively, similar effects can be reached 
by artificial inclusion of small but non-zero initial trans- 
verse flow. We demonstrate that such an imitation of 
the shear viscosity indeed improves the V2 coefficients 
and R ut/Rside ratios. 

The hadron yields and particle number ratios for most 
abundant hadrons measured by ALICE are also well de- 
scribed and explained in the model. The contributions of 
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Figure 7: (Color online) ir~ ,K~ and proton spectra and elliptic flow of all charged hadrons at midrapidity for 2.76 TeV LHC 
energy and different centrality classes (top) and HBT radii of 7r~ pairs for most central events, calculated in the "hybrid" 
model and compared to experimental data from ALICE [3,113, j4j|. Solid lines correspond to IC without transverse flow, while 
dashed lines stand for small initial transverse flows, (vt) =0.014 - 0.018. Corresponding HBT radii for the top RHIC energy 
are shown for comparison purposes for the basic case without initial transverse flow. Note that calculations for c = 10 — 20 % 
are not provided, but the ^-coefficient for c = — 5 % is presented instead as prediction. 
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different processes at the hadronic cascade stage to the 
final hadron yields are found. 

Note that the most of the parameters fixed from agree- 
ment with RHIC data are applied also for the description 
of 2.76 TeV Pb+Pb LHC collisions with different central- 
ity classes. In fact, the only changes are the normaliza- 
tion of initial entropy density distribution, contribution 
from binary collisions and baryonic chemical potentials 
at chemical freeze-out, which become zero. Then pion, 
kaon and (anti) proton spectra, p^-differential elliptic 
flow coefficients, and HBT radii are reproduced well also 
at LHC. It is the evidence of the similarity of the soft 
physics at the RHIC and LHC energies. 

As the next step it is planned to include some missing 
features of the hydrokinetics into future versions of the 
model, namely the prethermal dynamics and viscosity. 
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